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^ : Abstract 

, In Type III seesaw model, there are tree level flavor changing neutral currents (FCNC) 

I in the lepton sector, due to mixing of charged particles in the leptonic triplet introduced to 



realize seesaw mechanism, with the usual charged leptons. In this work we study these FCNC 
effects in a systematic way using available experimental data. Several FCNC processes have 



Oh; 

O I been studied before. The new processes considered in this work include: lepton flavor vio- 

lating processes r — > P/, r — > VI, V ^ W , P ^ W , M — > M'W and muonium-antimuonium 
■ oscillation. Results obtained are compared with previous results from li Ijhh, h Ijl^ 

(N : 

00 . Z ^ II and ji — e conversion. Our results show that the most stringent constraint on the 

o 

'nJ" ' e-to-r FCNC effect comes from r — > ir^e decay, r — > p'^/x and r ir^fi give very strin- 

' gent constraints on the /i-to-r FCNC effect, comparable with that obtained from r — > /ififj, 

o^ : 

, studied previously. The constraint on the e-to-fj, FCNC effect from processes considered in 

^ this work is much weaker than that obtained from processes studies previously, in particular 

X' 

^ • that from /i — e conversion in atomic nuclei. We find that in the canonical seesaw models 

the FCNC parameters, due to tiny neutrino masses, are all predicted to be much smaller 
than the constraints obtained here, making such models irrelevant. However, we also find 
that in certain special circumstances the tiny neutrino masses do not directly constrain the 
FCNC parameters. In these situations, the constraints from the FCNC studies can still play 
important roles. 
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Introduction 

Neutrino oscillation experiments involving neutrinos and antineutrinos coming from 
astrophysical and terrestrial sources have found compelling evidence that neutrinos 
have finite but small masses. To accommodate this observation, the minimal standard 
model (SM) must be extended. Some sensible ways to do this include: (a) Type I 
seesaw with three heavy right-handed (RH) Majorana neutrinos (b) the use of 
an electroweak Higgs triplet to directly provide the left-handed (LH) neutrinos with 
small Majorana masses (Type II seesaw |^), (c) introducing fermion triplets with zero 
hypercharge (Type III seesaw 3|), (d) the generation of three Dirac neutrinos through 
an exact parallel of the SM method of giving mass to charged fermions, and (e) the 
radiative generation of neutrino masses as per the Zee or Babu models ^. But in 
the absence of more experimental data, it is impossible to tell which, if any, of these 
is actually correct. Different models should be studied using available data or future 
ones. In this work, we carry out a systematic study of constraints on possible new 
flavor changing neutral currents (FCNC) in Type III seesaw model. 

The fermion triplet S in Type III seesaw model transforms under the SM gauge 
group SU{3)c X SU{2)l x U{1)y as (1,3,0). We will assume that there are three copies 



of such fermion triplets. The model has many interesting features 



5j, including the 



possibility of having low seesaw scale of order a TeV to realize leptogenesis iGl and 
detectable effects at LHC [3, 18|] due to the fact that the heavy triplet leptons have 
gauge interactions being non-trivial under the SU{2)l gauge group, and the possibility 
of having new tree level FCNC interactions in the lepton sector [ol, Q, H]. Some of 



the FCNC effects have been studied, such as U — > Ijlkk, U ~* ^jl: Z — > ijj and /i — e 
conversion processes. Several other FCNC processes studied experimentally have not 
been studied in the context of Type III seesaw model. We will study constraints on 
FCNC related to charged leptons in a systematic way using available experimental 
bounds listed in Ref. [12] by the particle data group. 

Before studying constraints, let us describe the model in more detail to identify new 
tree level FCNC in the lepton sector. The component fields of the righthanded triplet 
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S are, 



V 



'- -NyV2 ) ' 



and the renormalizable Lagrangian involving S is given by 



where Ll = (j^lJl)'^ is the lepton doublet. H 



kO\T 



+ h + iv)/V2f 



is the Higgs doublet with v being the vacuum expectation value, and H = iT2H*. 

Defining E = Ej^"^ + E]^ and removing the would-be Goldstone bosons t] and 0^, 
one obtains the Lagrangian 

£ = EipE + N^ipN'^ - EM^E - (n^—N^^ + h.c. 



+ 9 [W;IKi,PrE + W+iV^VL^ + h.c. ~gW;^E^,E 



-j={v + /i)iV^Fs^^L + {v + h)EY^lL + h.c. j . 



(3) 



One can easily identify the terms related to neutrino masses from the above. The 
mass matrix is the seesaw form 



Yj:^v/2V2\ / 



h.c. . 



(4) 



\^ Yj:v/2V2 Me/2 

The charged partners in the triplets mix with the SM charged leptons resulting in 
a mass matrix of the following form 



I, 



Yj:v Mj,) [El 



h.c. . 



(5) 



One can diagonalize the fermion mass matrices and find the transformation matrices 
between fields in weak interaction basis and in mass eigenstate basis defined as 



L,R 



u. 



( I' 



L,R 



( 



L,R 

R^ 



( 



U 



(6) 



3 



where the primed fields indicate mass eigenstates. Ul.r are (3 + 3)-by-(3 + 3) matrices 
if 3 triplets are present, and can be written as 

/ 



U, ^ I I . u, 

Ulei Ulee 



Unii Urie 
Urei Uree 







U.N 1 








\ Unu 


Unn I 



(7) 



To order v'^M^'^, one has 
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Uui = l-e, Ulie = Y^M^'v, Ulei = -M^'Yj^v , Ulee = I - 
Umi = 1 , Urie = miY^M^^v , Urei = -M^^Y^niiv , Ulee = 1 
U,, = (1 - e/2)UpMNS , U,N = Y^M^'v/V2 , 

Unu = -M^'Y^U,,v/V2 , Unn = 1- e'/2 , 



where 

e = F^Ms'1ei;V2 = U.nUI^ e' = M^^Yj^Y^M^^^ / 2 = uI^U.n ■ (9) 

Here Urmns denotes the lowest order Pontecorvo-Maki-Nakagawa-Sakata (PMNS) 
mixing matrix which is unitary. We have kept higher order corrections to the Up,, 
matrix. 

Using the above, one obtains the couplings of Z and physical Higgs h to the usual 
charged leptons 



Cnc = -^h^ {Pl{-\ + sin^ ew-e)+ Pr sin^ Ow ] IZ, 



cosOw V 2 



M 1 



'^H = TTTT-'^ (PLnii (3e - I) + Pr (3e - 1) nii) Ih . (10) 

ZIVi\Y 

Here we have dropped the "prime" on the fermion mass eigenstates. e is a 3-by-3 
matrix. Non-zero off diagonal elements in e are the new sources of tree level FCNC in 
charged lepton sector. The Z and Higgs coupling to quarks are the same as in the SM. 
We will use available FCNC data in a systematic way to constrain the parameter eui. 

Several processes, such as /j Ijlkh-, h ^1-, Z —* IV and fi — e conversion in 
atomic nuclei, have been studied and stringent constraints have been obtained for 
eiii which will be used as standards for constraints obtained from new lepton fiavor 
violating (LEV) processes considered here, r PI, r VI, V If, P II, 



M — M'lV and muonium-antimuonium oscillation. It turns out that with currently 
available experimental data, the LFV processes considered in this work involving r 
leptons provide very stringent constraints on the FCNC parameter eir- Our results 
show that the most stringent constraint on comes from r vr^e decay, r p^fi 
and r n^fi give very stringent constraints on e^^, comparable with that obtained 
from r fififi in previous studies. The strongest constraint on comes from fi — e 
conversion in atomic nuclei studied previously. We now present some details for the 
new processes mentioned above. 



Constraints from t ^ PI and t ^Vl 

Exchange of Z boson between quarks and leptons can induce t ^ PI and t ^ VI 
at tree level, where a pseudoscalar meson P = vr*^, r/, r/' and a vector meson V = p°, a;, 
and a charged lepton I = e, fi. The decay amplitudes for r Ml (where M denotes 
either V or P) can be written in the following form 

M = 2V2GFeir E {M{pJ\q-fM3PL-QqSm^ew)q\0)-\l{p,)r{l-l5)r{Pr 



q=u,d,! 



q=u,d,s 



where Gp is the Fermi constant, Qg is the electric charge of q-quark in unit of proton 
charge. J3 = 1/2 and —1/2 for up and down type of quarks, respectively. The factor 
9y = \~^ sin^ 9w and g'l = —\ for up type of quarks, and g'^, = ~| + | sin^ 9]v and 
g\ = \ for down type of quarks. The Pri Pi and p^^- are the momenta of r, / and M, 
respectively. 

For 7r°Z, the decay constant is defined as 

{n^{pn)\uiclM^) = -{'^°{Pn)\d'Jal5d\0) = -i-^{p^)o, (12) 

and its value is = 130.4 MeV. For r~ rjl and rj'l, due to the r] — t]' mixing, 

the decay constants /^(,) and /^(,) are defined as 

{v^'\Pnio)\s-falrM^) = -if^(o{P^('))a, (13) 
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where 



= ^sin^s + ^cos^o, /,^' = -2^sin^8 + ^cos^o, (14) 



with /s = 168 MeV, /o = 157 MeV, = -22.2°, and Oq = -9.1° (l3|. 



For r — > V/ decays, the decay constants /p, /^^ and are defined by 



(p°(Pp) |w7aw|0) = -{p\p,)\d^ad\0) 



V2 



m„ e 



{u{p^)\u-fau\0) = {uj{p^)\d-fad\0) = -^^^{eS)a, 

(0(p</.)|s7as|O) = Um^ie^)a, (15) 

where {ev)a is the polarization vector of V. We use fp = 205 MeV, = 195 MeV and 
/0 = 231 MeV Q. 

Exchange of Higgs boson can also induce qq couphng to II'. However, Higgs-mediated 
diagrams do not contribute to r — > P/ and t ^ VI because the bi-quark operator in this 
case is of the form qq which induces a vanishing matrix element for < P (or V^)|gq'|0 >. 

The decay rate for r~ ^ PI {P = 7r°, t], t]', and I = e, fi), averaged over the spin of 
r and summed over the spin of /, is given by 



ry2 f2 

^FJ Pi |2 I ^ I 
aPT^ 2 \^lr\ \Pi\ 

zirmi 



+ Till — 2m; — (m^ + m^)mp) , (16) 



where |pj = y {fnl + — mfy — Am'^m'p / {2mT-) . In the above expression, the decay 
constant fp is given by fp = f^^ with ap = 1 for r~ — > vr"/, and fp = /^(,) with 
ap = 1/2 for ^ rj W/. In the case of r — > rj^'H, the u and d quark contributions 
to the matrix element {i]^'^q'Ja'y5q\0) cancel each other in Eq. (fTTl) so that only the s 
quark contribution to the decay constant, /^(,), remains. 

Similarly, the decay rate for ^ VI {P = u, (p, and I = e, /i) is given by 



G pfyTny |2 I ^ 

O'V 2 — \Pi 



1 



2 I 2 2 I / 2 I 2 2\/ 2 2 2 

+ 777,; — my + —Y[rn^ + my — m^ j(m^ — — m^ 



(17) 



where |pj = y {jnl + my — m^)^ — 4m^my/(2m^). The decay constant fy is given by 
fy = fp with ay = (1/2 — sin^ OwY for r~ p^l, and /y = f^ with ay = (sin^ 6'vi//3)^ 
for ^ u;/, and fy = f^j, with ay = 2(1/4 — sin^6'^y/3)^ for r~ — 0/. 
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TABLE I: Constraints from t ^ PI and t ^Vl. 



Process 


Branching Ratio 


Constraint on e 








< 8.0 X 10-^ 


-8 


1 ^eT 1 


< 4.2 X 10- 


-4 


t" 




< 1.1 X 10-^ 


-7 


1 ^M''' 


< 7.0 X 10' 


-^A 


T 


— > r]e~ 


< 9.2 X 10" 


-8 


1 ^eT 1 


< 1.2 X 10" 


-3 


T 


^ W~ 


< 6.5 X 10- 


^8 


1 ^AtT 


< 9.7 X 10' 


-4 




^ rj'e~ 


< 1.6 X 10~ 


^7 


1 ^eT 1 


< 1.0 X 10- 


-3 


T" 




< 1.3 X 10-^ 


-7 


1 ^/XT 


< 1.0 X 10- 


-3 






< 6.3 X 10" 


-8 


1 ^eT 1 


< 6.5 X 10- 


-4 


T~ 


- - pV- 


< 6.8 X 10- 


^8 


1 ^/XT 


< 6.8 X 10" 


-4 


T 


— > tje^ 


< 1.1 X 10- 


-7 


1 ^eT 1 


< 3.2 X 10- 


-3 






< 8.9 X 10- 


-8 


1 ^/XT 


< 2.5 X 10" 


-3 


r 


— > (j)e~ 


< 7.3 X 10" 


'8 


1 ^eT 1 


< 7.5 X 10- 


-4 




- (p^i- 


< 1.3 X 10- 


^7 


1 ^/XT 


< 1.0 X 10- 


-3 



Using the current experimental bounds on the branching ratios, we find the 
constraints on the parameters legrl and |e^T-| which are shown in Table [B Notice that 
the constraint on |eer| obtained from 7r°e~ is jeerl < 4.2 x 10~^, which is more 

stringent than the so far most stringent bound obtained from r eee as shown in 
Table IIVI The constraints on \e^r\ obtained from r~ ir^fi^ and t~ p^^~ are 
comparable to the so far most stringent bound shown in Table IIV[ The upper bounds 
on |ee(/x)T| from r~ rj^'H and r~ ujI are weaker. 

Constraints from V II' and P II' 

Here V can be a vector meson J/ip or T, and P can be a pseudoscalar meson vr*^, 77 or 
rj'. The / and /' stand for charged leptons with different fiavors / 7^ These processes 
can be induced by exchange Z boson between quarks and leptons. The general decay 
amplitude for M — II' (where M denotes either V or P) is given by 

M = 2V2GF ew {^IiiMsPl - Q,sm^ew)q\M{pJ) ■ [ui{pi)ni - 75)t'r(p2)] 

q=u,d,s,c,b 
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1^p2Gp eiv 

q=u,d,s,c,l 



where we use the decay constants /j/^ = 416 MeV and fr{3S) = 430 MeV [1^, Il5 |. 
Again, exchange of Higgs boson does not contribute to these two classes of processes 
since (0|gg|M) = 0. 

The decay rate for V — > iT {V = J , T) is found to be 



2 1 1 



1 



[nil - mi,) 



2ml 



(19) 



'my + mf — mf,y—4mYm'i/{2mv), and g'^ = g'^ for V = J/\E' and 



where |pj 

Similarly the rate of a pseudoscalar meson decay P II' {P = tt^, t], t]') is given by 



'1 = for V = T. 



ry2 f2 

■p ^FJP\ \2 1^1 [/ 2 I 2\ 2 / 2 2\2 

T = ap- \eu'\ bj [mi + mi,)mp - [mi - mi,) 

znm p L 



(20) 



where 



m. 



ml 



— mf,)"^ — Ampmf /{2mp), and ap = 1, fp = for P 



and ap = 1/2, fp = /^(,) for P = r]^'\ Note that as in the case of r~ rj^'^l, only the s 
quark contribution to the decay constant, /^(,), appears in r]^'^ I I'. We find that the 
constraints on \eii, \ from these two body meson decays are rather weak as summarized 
in Table [Tll The constraints obtained are much weaker than those obtained in the 
previous section. 



Constraints from M M'W 

We now consider semileptonic three body decays of the type M 



M'W with 



M = B,K and M' = K, K* tt, such as B ^ Kll', B K*ll', B 



TT 



//', and K nil'. 



These decays can occur through quark level subprocesses h — * slV 01 s dlV . The 
FCNC h ^ s 01 s ^ d transition can arise via Z-penguin and Higgs-penguin diagrams 



at one loop level the same way as in the SM. After taking into account t 
h-s-Z and 6-s-Higgs couplings (or s-d-Z and s-(i-Higgs couplings) [la. 



le SM effective 



17|, the lepton 



flavor violating FCNC processes h — > slV (or s dlV) can occur at tree level via the 
couplings given in Eq. (fTOl) . 
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TABLE II: Constraints from V IV and P ^ IV . 



Process 


Branching Ratio 


Constraint on e^^, 


T(35) - 




< 5 X lO"'^ 


1 ^er 


< 0.39 


T(35) - 


^l^T^ 


< 4.1 X 10-6 


1 ^^J-T 


< 0.35 


J/*(15) 




< 1.1 X 10-6 




~ 0(1) 


J/*(15) 


^ e^T^ 


< 8.3 X 10-6 


1 ^ET 1 


~ 0(1) 


J/^{IS) 




< 2.0 X 10-6 


1 ^IJ'T 1 


~ 0(1) 






< 3.4 X 10-9 


1 ^e/x 


< 0.80 






< 3.8 X 10-1° 


1 ^e/^ 


< 0.27 






< 6 X 10-6 




~ 0(1) 






< 4.7 X 10-^ 




~ 0(1) 



The decay amplitude for M M'lV [where M = B {or K); M' = K,K*,7i (or vr); 
/, /' = e, /i, T {I ^ I'))] is given by 



(21) 



where and Ai^ denote the Z-mediated and Higgs-mediated decay amphtude, re- 
spectively, in the following form 



:V*„,iVini- 



327r2 cos^ OwM^ 

X [m/(A:i)7"(1-75H'(A:2)] 



Co{xi) eu> (M'(p')|g"7a(l-75)glM(p)) 



(22) 



9 _ ^rnlmq: 



eu> {M'{p')\q"{l+^,)q'\M{p)) 



10247r' 

X {ui{kx)[{mi + mv) + (m^ - mi)'^^]vv{k2)} 



(23) 



where (i) for B I&'^IV , q' = b and q" = s, (ii) for B nil', q' = b and q" = d, 
(iii) for K nil', q' = s and q" = d. The Viqi denotes the CKM matrix element 
with i = t,c,u and Co(xi) = (xj/S) [{xi — 6)/(xi — 1) + (Sxj + 2) lnxj/(xj — 1)^] with 
Xi = mj/mly |l6|. 

Compared with the Z-mediated amplitude, the Higgs-mediated amplitude is negli- 
gibly small, since rrih ^ m^, mi, so that the Higgs contribution can be safely neglected. 
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For example, in the cases of -B — > K^*^lV and K — > Till' decays, M^'' j M.^ is suppressed 
roughly by 0{xt{m{,mi/mD) and 0{xt{msmi/mD), respectively. 

For B Pll' (P = IT, K), the form factors Fi and Fq (or /+ and /_) are defined by 



(P(p')|^7.(l-75)6|S(P)) = F,iq') 



{P + P')a - 



— m 



K 



— m 



K 



(24) 



where q = p — p'. For B K*ll', the form factors V, Aq, Ai, and A2 are defined by 

2V{q') 



{K*{p',e)\s^^{l-^,)b\B{p)) = -e^p^^e^^pfp'-- 



-i 



■ q 



rriB + rriK* 
Qa {rriB + 'mK')Ai{q'^) 



+i y{p + p)a - 

2mK*{e* -q) 



-Qa (e* ■ q)- 



—i- 



niB + rriK* 

(25) 



where e is the polarization vector of the K* meson. For numerical anal ysis , we use 
the form factors calculated in the framework of light-cone QCD sum rules 18||. The q^ 
dependence of the form factors can be expressed as 

m 



F{q' 



2 ' 



(26) 



where the values of the parameters F{0), ap and bp ioi B it, B ^ K and B — » K* 
are given in fl8|. 

Summing over the spins of the final leptons, we obtain 
dT{B Pll') 1 



Cl{x,)\eiv 



dq' 



1927r5 sin^ Ow cos^ Ow 



A^/^(m|, rrip, q^^ 



m 



B 



1 + /+(?') +3p /o(g 



x(l-2p)2 

where A(a, 6, c) = (a — 6 — c)^ — 46c, p = m//(2g^) and 



(27) 



(28) 



{m%, mp, q^ 

Here the mass of only one light lepton in the final state has been neglected so that 
the parameter p represents the effect of the remaining lepton mass, e.g. rrir- Thus, for 
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B Kefi decays, p can be neglected. The decay rate for B K*ll', summed over 
the spins of the final leptons and K*, is given by 

dT{B K*ll') 1 



ds 



7687r^ sin'' 6w cos^ 9w 
8m%sX 



X 



{tub + mK*y 



+ 



12s 



-2mB Re 



4 

TUB 



(niB + rriK* 



A(l - r - s) 



(29) 



where r = m|-,/m^, s = g^/m^, and A = 2r— 2s— 2rs. The branching ratios 

for B — s> ttU', B — s> Kll' and B K*ll' can be calculated after the decay rates given in 
Eqs. (1271) and fl29|) are integrated in the range {mi + m//)^ < Q'^ < {'^b ~^m'Y- From 
the current experimental bounds on those branching ratios, we obtain the constraints 
on shown in Table IIIII 

For K — s> Txll\ we normalize the branching ratio to Tr^e^Ve and neglect the 

phase factor difference [19|]. We have 



B{K+ ^ 7r+/P) 2a2 V^^Vm ' 



B{K+ - 
B{Kl 



vr^ sin 9\y cos"^ 9w 



Vu. 



Cl{xt) \ew 



2a' 



B{K^ 



tk+ vt^ sin 9w cos^ 9 



w 



Im 



(30) 



where tk is the lifetime of the Kaon. Note that the mo del- dependent form factors do 



le above formulas. Using the experimental value B{K^ 



12| , we obtain the constraints on ew shown in Table IIIII Alternatively, 



not appear in t 
(5.08 ±0.05)% 

the decay rate for K nil' can be calculated by using Eq. (I27|) . In this case, the mass 
of muon is not neglected and the parameter p = m^/(2g^). The relevant form factors 
are given by 



-1 - \W 



f^\q') ^ f^^q') + 



m 



K 



m^2 



f^^q') ^ -1 - Aog^ 



(31) 
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TABLE III: Constraints from M M'lV . 



Process 


Branching Ratio 


Constraint on e 


i' 


B+ 




< 6.4 X 10' 


-3 


1 ^e/i 1 


~ 0(1) 




B+ 


— > vr+e"//"*" 


< 6.4 X 10" 


-3 


1 ^e/^ 1 


~ 0(1) 




B+ 


— > 7r"^e^/i^ 


< 1.7 X 10- 


-7 




< 0.56 




B+ 




< 9.1 X 10" 


-8 




< 0.18 




B+ 




< 1.3 X 10" 


-7 




< 0.21 




B+ 




< 9.1 X 10' 


~8 




< 0.12 




B+ 




< 7.7 X 10' 


-5 


1 ^f^T 1 


~ 0(1) 








< 1.4 X 10' 


-7 




< 0.73 




5° 




< 2.7 X 10' 


-7 




< 0.21 




B+ ^ 


K* (892)+ e+fi- 


< 1.3 X 10" 


-6 


\eef,\ < 7.1 X 10- 


-2 


B+ ^ 


A'*(892)+e^/i+ 


< 9.9 X 10' 


-7 


\eei,\ < 6.2 X 10' 


-2 


5+ ^ 


i^*(892)+e±/i^ 


< 1.4 X 10' 


-7 


\eet,\ < 1.7 X 10' 


-2 




i^* (892)0 e+/i- 


< 5.3 X 10' 


'7 


|ee^| < 4.5 X 10' 


-2 


B^^ 


i^* (892)0 e^/i+ 


< 3.4 X 10' 


-7 


\eef^\ < 3.6 X 10' 


-2 




A'*(892)0e±/i^ 


< 5.8 X 10' 


-7 


\eef,\ < 3.4 X 10' 


-2 


K+ 




< 1.3 X IQ- 


11 


lee^l < 0.44 [0.^ 




K+ 


— > 7r"'"e~/i~'" 


< 5.2 X 10" 


10 


1 ^e/^ 1 


~ 0(1) 




Kl 


^ 7rOe±;uT 


< 6.2 X 10' 


-9 


1 ^e/^ 1 


~ 0(1) 





y. 



where A+ = 0.067 fm^ and Aq = 0.040 fm^ |20|. The constraints on ew obtained in this 
way (number shown in the bracket for 7T~^e~^fi~) is similar to those obtained by 

using Eq. ( 1301) as shown in Table lllli The constraints obtained here are again much 
weaker than those obtained from t ^ PI. 



Constraint from muonium-antimuonium oscillation 

At tree level, exchange of Z boson can generate an effective Hamiltonian of the form 

H^jf = V2GFe:lfij,{l - 75)e/x7'^(l - 75)e . (32) 
12 



TABLE IV: Constraints from /j — > Ijlkk, h — > hi decays and /i — e conversion. 



Process 


Conversion rate 


Constraint on e 


I' 1 




- e 


conversion 


< 4.3 X 10-12 




< 1.7 X 10" 


-7 


Process 


Branching Ratio 


Constraint on e 


i ' 






-> e" 


e 


< 1 X 10-12 


1 ^e/x 1 


< 1.1 X 10" 


-6 


T~ 




e" 




< 3.6 X 10-8 


1 ^ET 1 


< 5.1 X 10" 


-4 


T~ 








< 3.2 X 10-8 


1 ^I^T 


< 4.9 X 10" 


-4 


T~ 








< 4.1 X 10-8 


1 ^ET 1 


< 7.2 X 10" 


-4 


T~ 








< 2.7 X 10-8 


1 ^fiT 


< 5.6 X 10" 


-4 








67 


< 1 X 10-15 


1 ^e/^ 1 


< 1.1 X 10" 


-4 








67 


< 5 X 10-11 


1 ^ET 1 


< 2.4 X 10" 


-2 










< 4 X 10-11 


1 ^A"" 


< 1.5 X 10" 


-2 



TABLE V: Constraints on e/// from Z IV decays. 



Process 


Branching Ratio 


Constraint on e^^, 


Z e±//T 
Z ^ e±rT 


< 1.7 X 10-6 

< 9.8 X 10-6 

< 1.2 X 10-5 


lee^l < 1-8 X 10-3 
\eer\ < 4.3 X 10-3 
\ei,r\ < 4.8 X 10-3 



This interaction will result in muonium-antimuonium oscillation. 

The SM prediction for muonium and antimuonium oscillation is extremely small. 
Observation of this oscillation at a substantially larger rate will be an indication of new 
physics. Experimentally, no oscillation has been observed. The current upper limit for 
the probability of spontaneous muonium to antimuonium conversion was established 
at Pmm < 8.3 X 10-11 (90% C.L.) in 0.1 T magnetic field [211 . 

In the absence of external electromagnetic fields, the probability Pmm of observing a 
transition can be written as [22] Pmm 

(OT) ~ |(5|V(2r2), where 6 = 2{M\Heff\M) and 
is the muon decay width. For H^jf given above, the transition amplitude is given by 
S = 32GFelJ{V27ia^) for both triplet and singlet muonium states, where a ~ {arrif,) i 
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is the Bohr radius. The probabihty Pmm has strong magnetic field dependence which 
usually occurs in experimental situation. With an external magnetic field, there is a 
reduction factor Sb, i-e. Pm]^i{B) = SbPmai{^'^)- The magnetic field correction factor 
Sb describes the suppression of the conversion in the external magnetic field due to the 
removal of degeneracy between corresponding levels in M and M. One has Sb = 0.35 
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231]. Using this experimental information, one obtains a 



for our case at i? = O.IT 
constraint 



lee^l < 4 X 10-^. (33) 

This constraint is rather weak compared with that from fi — e conversion. 

Exchange of Higgs boson will also contribute. But this contribution is suppressed 
by a factor m?^/m\ and can be safely neglected compared with Z boson exchange 
contribution. 



Constraints from /j — > Ijlkhi h ~^ Z II' decays and M — e conversion 

These processes have been studied in the literature before [o], [l^. For comparison, 
we summarize the results for constraints on eur for Ijlkk, h — Ijl and /i — e 

conversion, and Z — > //' 2J] in Tables IIVI and |Vl respectively. The most stringent 
upper bound on |ee^| is of order 10"'' from yU — e conversion in atomic nuclei. The upper 
bounds on |eeT-| and |e^^| obtained are of order 10~^ from r — s> eee and r — > /i/i/i decays. 



Discussions on the mixing matrix U^n between the light and heavy neutrinos 

We now discuss some implications of the constraints obtained earlier on the model 
parameters. In this model, to the order we are studying, the light neutrino mass is 
related to U^n with 

t^PMNS'^i/f^PMNS = —Uyf^M^JJl^ , (34) 
where the light neutrino mass matrix rhy is diagonal: 

rriy = diag (m^,^, m^^, m^g) = f/pMNs"^i.f^PMNS • (35) 
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Thus, one might think that the elements of Uun are too small to be relevant to the 
FCNC discussion, because with only one generation of the light and heavy neutrinos, 
\Uun\ is simply given by (m^/Ms)-'^/^. It leads to the fact that for Ms > 100 GeV, Uun 
is less than 10~^, since the light neutrino masses must be less than an eV or so. If with 
more than one generation of the light and heavy neutrinos, all elements of U^n are the 
same order of magnitudes (the canonical seesaw models), the resulting elements of the 
e matrix will all be way below the constraints we have obtained. This makes the model 
irrelevant for an experimental detection. The FCNC study of the kind studied here is 
therefore not interesting for canonical seesaw models. However, it has been shown that 
with more than one generation of the light and heavy neutrinos, there are non-trivial 
solutions of f/j^Tv such that the right hand side of Eq. flMl) becomes exactly zero but 



the elements of U^n can be arbitrarily large 25|, |26|] . Thus, these solutions evade the 
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canonical seesaw constraint \U^n\ = {ttT'u/Mj^Y^'^ held in the one generation case 
261 ]. It is interesting if one can find the Uun which satisfies existing experimental 
constraints by adding small perturbations to the above non-trivial solutions. A recent 
study has shown such solutions of U,yN that indeed can have large elements and satisfy 



the current experimental constraints [26|. In the following we will describe some of 
those solutions having relevance to our FCNC study. 

Let us indicate the solution of 11,,^ which gives zero light neutrino mass as Uq. We 
then add a perturbation Us to Uq such that U^n = Uq + Us- Since UqMj^Uq = 0, the 
neutrino mass matrix is given by 

= -f/oMsf/J - UsM^U^ - U^M^Uj . (36) 

If the first two terms are not zero, the matrix elements aij in Uq and 6ij in Ug are of 

order aij6ij ~ rriu/Mj] which is much smaller than 1. Since we are interested in having 

large atj, the elements 6ij must be much smaller than aij, and the third term, for 

practical purpose, can be neglected. If on the other hand, the first two terms are zero, 

the third term must be kept. The elements of Us in this case are of order (mj^/Ms)^/^. 

In the basis where Ms is diagonal, one can write 

/ 1 1 1 \ rriM , , 

M. = Ms = diag(^-,-,-jm,, = — ^ (37) 
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where, for convenience, we have introduced a scale parameter ttt-at to represent the scale 
of the heavy neutrino, which we choose to be the lightest of the heavy neutrinos. The 
contribution to e is given by 



UunUIn 



(3J 



We show three types of solutions relevant to our study of FCNC: (a) sizeable 612,13,23; 
(b) sizeable 623 and small 612,13; ^i-iid (c) sizeable 613 and small 612,23- In case (a), the 
data from fi — e conversion in atomic nuclei constrain |6i2| to be less than 1.7 x 10~^ 
which makes 613,23 too small to be of interest. We therefore need to find other classes of 
solutions where 612 is automatically much smaller than 613,23. These are the cases (b) 
and (c). If these types of solutions are correct, the constraints from r decays discussed 
previously in this paper are still relevant for experimental search. 



The numerical results will be given by using the central values of Am 
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(7.65i[];^g) X 10-5 eV2 and jAm^^l = {2A0+l\j) x IQ-^eV^, determined by a recent fit 



to global neutrino data |27|], and t/pMNS in the tri-bimaximal form 28| for simplicity 

\ 



tribi 



1 1 1 

Ve Vs V2 

\ ^ ^ ^ I 

\ Ve Vs V2 / 



For the details of the following solutions, we refer to Ref. 
For desired solution is given by 

^ a a i\pla ^ 
f^PMNS h b i\plh , = ?7pmns 
y c c i\/2c I 



(39) 



26| 



•^0 



I (5ii (5i2 (5i3 ^ 

<^21 <^22 

\ h\ h2 ) 



n 



(40) 



where TZ = diag V^j ■ There are two types of solutions corresponding to 

normal and inverted hierarchies in light neutrino masses, but always one of the masses 
becomes zero as follows, 
(i) Normal hierarchy: 

a = 0, rhj, = diag ( 0, — 1, — 1 256mAr , 
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0.00050 




FIG. 1: For case (a), the upper limits on the magnitude of the element of UyN in terms of 
the heavy neutrino mass parameter r = ri + r2 + 2r^ . The solid and dashed lines correspond 
to the normal and inverted hierarchy cases, respectively. 



0.33 0.33- 0.97i 0.33 + 0.97i^ 
0.33 + 0.97i 3.1 -2.5 + 1.9i 

0.33 - 0.97i -2.5-1.9i 3.1 j 



|6|'r 



(41) 



(ii) Inverted hierarchy: 



m,. 



0.99 



diag f — , -1, j 256m AT 



0.01-0.70i 0.01-0.70i 



0.01 + 0.70 i 
0.01 + 0.70i 



0.50 
0.50 



0.50 
0.50 



|6p r 



(42) 



where 5 = 52i+S22+iV^523 and r = ri+r2+2r3. From /i — e conversion in atomic nuclei 
(1^12! = l^e/il < 1.7 X 10"''), \b\^/r is constrained to be smaller than 4.1 x 10~^ (normal 
hierarchy) or 4.9 x 10"'^ (inverted hierarchy). In both cases, 1 613 23] are constrained to 
be less than 0(10^'') which are way below the best constrained from r fiflfi and 
r — >• 7r°e decays. 

In Fig. [1] we show the upper limits from the /i — e conversion constraint on the 
magnitude of the element b of U,^n in terms of the heavy neutrino mass parameter 
r. Since is the lightest of Mi, r is in the range 1 < r < 4. Depending on 
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the heavy neutrino mass hierarchy, the value of |6| can be different. With the same 
constraint, to have the largest b, one would require the two heavier ones to be much 
larger than the lightest rriN- As far as FCNC processes are concerned, the hierarchy of 
the heavy neutrinos is not important because the parameter always involves r. But for 
the production of a heavy lepton at LHC, via qq' — ^ W* IN or qq — > {Z*, h*) ^ IE 
for example, it is preferred to have a larger 6, because in that case, not the combination 
+ ^2 + 2r3) but the individual \h\'^rj is relevant to the production cross section. 



For case (b), the following form serves the purpose with the choice 



uN 



Uip^ + Ul where 







7^ 



a ia 
\ h ih ) 

Here a is of order [(a, h)5i 



a 





\^ y 



7^ 









Vo o; 



7^ . (43) 



11/2 



so that one should keep terms in the calculation. 



neglecting 5ij5ki and a5ij terms. The eigen-masses are 



diag(a 5 



12 



a' 



-2a — o? ■, 0)m^ 



(44) 



and so this is an inverted hierarchy case with rny^ 
given by 

/^O 0^ 



0. Numerically, the matrix e is 



1 1 
1 1 



where p 



r^. Thus, the constraint |e 



23 1 



< 4.9 X 10~^ from r 



(45) 



/i/i/i 



decays translates into |a|^/p < 2.2 x 10~^. Since ri does not show up in Uq in this 
case, it would be more convenient to choose ttt-at to be the lightest of M2,3. 

For case (c), the desired results can be obtained by choosing U^,^ = Uq+ U^^^ + [/^, 



with 



a ia 

b ib 



n 



a 
/? 




n . 



(46) 



0.022 



0.020 



Case (b) 



0.018 



Case (c) 



0.016 - 



1.0 



1.2 



1.4 1.6 
P 



1.8 



2.0 



FIG. 2: For cases (b) and (c), the upper limits on the magnitude of the element of U^n in 
terms of the heavy neutrino mass parameter p = r2 + r^. The sohd hne corresponds to case 
(b) and the dot-dashed and dotted hnes correspond to the normal and inverted hierarchy 
cases, respectively, in case (c). 

This particular choice allows all the three light-neutrinos to have nonzero masses. 
Taking = 0.1 eV, two possible solutions are found and give the matrix e as follows. 

(i) Normal hierarchy (with m^j^ = 0.0996eV and m^.^ = 0.111 eV): 

'^ 1 0.001 - 1.0i^ 



0.001 + 1.0z 1.1 

(ii) Inverted hierarchy (with rriui = 0.0996 eV and mi,^ = 0.0867 eV): 

^ 1 0.001 + 0.96i^ 



0.001 - 0.96i 0.93 



V 



(47) 



V 



(48) 



The bound |e 



131 



< 4.2 X 10"^ from r 



TT 



'e decays then implies |a|i/p < 



2.0 X 10~^ in the two cases. 

In Fig. [2] we display the upper limits on the magnitude of the element a of Uun in 
terms of the heavy neutrino mass parameter p for cases (b) and (c). In this case, p 
is in the range 1 < p < 2. With the same constraint, the hierarchy that the heavier 
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of M2,3 is much larger than tjin would be required to obtain the largest a. Similarly 
to case (a), concerning FCNC processes, the hierarchy of the heavy neutrinos is not 
important. But, concerning the production of a heavy lepton or at LHC, a large 
cross section can be obtained for ^ 115 GeV [26 1. 

The above examples clearly show that with the constraints from FCNC transitions 
as well as from the tiny neutrino masses, the elements of 1/^^ can still be large. There 
is another class of processes which also provides constraints on the elements of U^n- 
These processes involve neutral currents conserving lepton flavor and can be used to test 
deviations from the SM predictions for electroweak precision data (EWPD) j^]- They 
have been measured mainly at LEP and provide bounds on the combinations of the 
diagonal elements of U^n. The constraints extracted from the EWPD are \{Ui^N)ii\ < 
(9(0.01) |29|. In contrast, the FCNC constraints discussed above involve combinations 
containing the off-diagonal elements and impose more stringent constraints, such as 
1^121 = \J2k{UuN)ik{U*^)2k\ < 1-7 X 10^'^. The non-zero elements of Uun in the two 
examples we give above with suppressed eu, being at most of (9(0.01), satisfy all these 
constraints. 

Large elements of U,^n also have important implications for a direct test of the model 
by producing the heavy neutrinos at LHC. The elements of U^n with the magnitude of 
order 0.01 are large enough to be detectable at LHC 26|. The heavy neutrino can 
be produced through the mixing via qq' — W* l^N. Similarly, the heavy charged 
lepton E can also be produced through the mixing via qq {Z*,h*) l^E^ and 
qq' — s> W* — > uE^. At LHC the production cross section for a single heavy neutrino 
A^ can be larger than 1 fb if the heavy neutrino mass is less than 115 GeV with the 
elements of U^n being 0.01. The production cross section of a single E is slightly 
smaller. This can provide useful information about this model. 



Conclusions 

We have systematically studied various FCNC processes in the lepton sector in the 
framework of Type HI seesaw model. Using the current experimental results, we have 
put the constraints on the parameters ew which are responsible for tree level FCNC in 
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the charged lepton sector. The new processes that have been considered are: the LFV 
processes r PI, r — > VI, V II', P II', M — * M'W and muonium-antimuonium 
oscillation. 

Although exchange both Z and Higgs bosons at tree level can induce FCNC in 
charged lepton sector, we find that there is no contribution from Higgs exchange in the 
processes r P{V)l and V{P) II', and the effects of Higgs exchange are negligibly 
small in the last two classes of processes. 

We now compare constraints on various FCNC parameters obtained from processes 
considered in this work with those obtained in previous studies. It turns out that 
with currently available experimental data, the LFV processes considered in this work 
involving r leptons provide very stringent constraints on the FCNC parameter e^^. Our 
results show that the most stringent constraint on comes from r — * vr^e decay with 
|eeT-| < 4.2 X 10~^. T p^fi and r ir^fi give very stringent constraints on e^^- with 
|e/ir| < 6.8x 10~^ and |e^^| < 7.0 x 10~^, respectively, comparable with |e^T-| < 4.9x10"^ 
obtained from r fififi in previous studies. The strongest constraint on comes 
from fi — e conversion in atomic nuclei studied previously with \eef^\ < 1.7 x 10^'. The 
new constraint on obtained from processes considered in this work is much weaker. 

Two body meson decays, such as T(35') ^ II' , J/^ ^ II', n ^ II' and 7]^''^ II', 
provide rather weak bounds on |e//'| at most of order lO"'^. The constraints from 
semileptonic three body B or K decays of the type M — > M'll' are also rather weak 
with upper bounds on \eii'\ in the range O(10~^) ~ 0(1). 

In the canonical seesaw models, where the elements of 11,^^ are of the same order of 
magnitude as that for an one generation seesaw model, {m,^/mj\iY^'^, it is not possible 
to have elements of e which are sufficiently large to reach the FCNC bounds studied 
in this paper. The FCNC effects studied are therefore not interesting for the canonical 
seesaw models. However, with more than one generation of light and heavy neutrinos, 
in certain special circumstances the mixing is not constrained directly by the tiny 
neutrino masses and therefore can be large. Thus in this class of seesaw models, it is 
possible to have large FCNC interactions. These circumstances have been studied by 



several groups 



25 



26|. We find some example solutions which can lead to the FCNC 



21 



parameters eij large enough to reach the constraints obtained here. The search for 
FCNC effects can still provide further information on the seesaw models. We comment 
that efforts in constructing models with certain symmetries to evade the canonical 
seesaw constraints on the mixing matrix UiyN have been made in various ways 30|]. It 
would be interesting to further study related phenomenology to test these models. 

We would like to comment that in some processes considered in this work it is 
possible to have CP violating signatures, such as lepton and anti-lepton decay rate 
asymmetries, and asymmetries in Z decays into //' and /'/ Sjj. To have non-zero 
effects, one needs not only a weak phase appearing in CP violating couplings coming 
from the complex e^j and Upmns matrix, but also a strong phase appearing in an 
absorptive part from loop induced decay amplitudes. Since we consider that the heavy 
seesaw scale M is heavier than Z, no absorptive part will be developed with the heavy 
triplets in the loop. Only light degrees of freedom in the loop for Z decays into W 
and /'/ can generate the absorptive parts which are generally small. The resulting CP 
violating effect will therefore be small. If polarizations of the initial and final particles 
can be measured, it is possible to construct CP violating observables which does not 
need the absorptive parts [32|. We will carry out the detailed studies elsewhere. 

Finally let us comment on several possible improvements on ew from future 
experiments. Improved data for t ^ PI and r ^Vl decays at various facilities, such 
as B and r-Charm factories, can improve the bounds on ew. Bounds from V W 
and P II' can also be improved, but may not be able to compete with constraints 
from other experiments. The current bound from B Kfir is rather weak. But 
at LHCb about 10^^ bb pairs are expected to be produced each year, and this decay 
mode may be useful in improving bound on e^^. Rare kaon decays will be studied at 
J-PARC with high precisions so that the current weak bounds from kaon decays may 
also become much stronger. But bounds obtained may still not be competitive with 
others, fi — e conversion in atomic nuclei will also be studied at J-PARC with several 
orders of magnitude improvement in sensitivity. Constraint on can be improved by 
more than an order of magnitude. It may be very difficult to improve constraint on 
from muonium-antimuonium oscillation to the level /i — e conversion can achieve. At 
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present Z II' do not provide the best bounds on ew. However, the Giga-Z modes at 
future colhders, such as ILC, the sensitivity can be improved by up to three orders of 



magnitudes 



331 ] ■ Future studies of Z — * er and Z fir may improve the bounds on 



eer and e^r- It is clear that FCNC effects in Type III seesaw model can be further tested. 
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